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Abstract 

The necessary and sufficient conditions for the exactness of the semiclassical approxi¬ 
mation for the solution of the Schrodinger and Klein-Gordon equations are obtained. It 
is shown that the existence of an exact semiclassical solution of the Schrodinger equation 
determines both the semiclassical wave function and the interaction potential uniquely up 
to the choice of the boundary conditions. This result also holds for the Klein-Gordon equa¬ 
tion. Its implications for the solution of the Wheeler-DeWitt equation for the FRW scalar 
field minisuperspace models are discussed. In particular, exact semiclassical solutions of the 
Wheeler-DeWitt equation for the case of massless scalar field and exponential matter poten¬ 
tials are constructed. The existence of exact semiclassical solutions for polynomial matter 
potentials of the form Xcj)^^ is also analyzed. It is shown that for p = 1, 2 and 3, right-going 
semiclassical solutions do not exist. A generalized semiclassical perturbation expansion is 
also developed which is quite different from the traditional h and My^-expansions. 
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1 Introduction 


The semiclassical or WKB approximation is usually discussed in textbooks on nonrelativistic 
quantum mechanics in the context of stationary states, i.e., determination of the energy eigenvalues 
and eigenfunctions, |^. This approximation can also be used to obtain approximate and in some 
cases exact solutions of the dynamical problem, i.e., full Schrodinger equation. To the best of my 
knowledge, however, the utility of the semiclassical approximation in obtaining exact solutions of 
the Schrodinger equation has not been fully explored. 

The same seems to be the case for the relativistic quantum mechanics. The importance of 
the semiclassical approximation in the relativistic case is probably best appreciated in quantum 
cosmology, , specihcally, in the analysis of the Wheeler-DeWitt equation which is essentially 
a Klein-Gordon equation on a superspace, @]. 

In more general terms, the semiclassical approximation is usually viewed as an approximation 
scheme in which one neglects all but the hrst term in an asymptotic perturbation expansion of 
the solution of a linear differential equation. Typical examples of such an asymptotic expansion 
are the loop expansions of quantum mechanics and held theory where the perturbation parameter 
is the Planck constant^] h, []^. In the context of quantum cosmology the relevant perturbation 
parameter is the gravitational coupling constant (or inverse of the Planck mass Mp), 

Usually, these perturbation expansions are singular and it is difficult, if not impossible, to obtain 
their precise structure. 

There is a more universal alternative for dehning the semiclassical approximation where the 
validity of the approximation is not linked with the values of the physical constants but determined 
by the properties of the wave function. In this approach one uses the polar representation of the 
wave function 

= ( 1 ) 

and obtains two coupled nonlinear differential equations for the amplitude R and the phase (angle) 

S' of "0 by substituting (|I]) in the dynamical equation. As it is demonstrated for the Schrodinger and 

^Here I have assumed that the kinetic term in the Lagrangian does not involve a coupling constant. For example 
in nonrelativistic quantum mechanics in a Euclidean space, this corresponds to the case where the mass m of the 
particle is set to unity. Otherwise, the perturbation parameter for the loop expansion is 
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Klein-Gordon equations in sections 2 and 3 below, there emerges a quantity called the quantum 
potential Q which controls the coupling of these two equations. In other words, if Q which 
depends only on R happens to be negligible, then one of the equations decouples from the other. 
The decoupled equation which only involves S turns out to satisfy a Hamilton-Jacobi equation. 
Thus, for Q = 0, S can be identihed with the classical action function of the corresponding 
classical theory. This observation is originally due to Bohm [Q. It provides the basis for the 
de Broglie-Bohm causal or ontological interpretation of quantum mechanics, [^]. The latter has 
recently been applied to problems of quantum cosmology by several authors, [0. 

The idea of the quantum potential leads to a precise criterion for the validity of the semiclassical 
approximation, namely the condition Q 0. More precisely, one has the following 


Definition: A wave function is said to he semiclassical if the corresponding quantum po¬ 
tential vanishes identically. 


Note that the quantum potential Q is determined by the amplitude R of the wave function. Thus, 
the validity of the semiclassical approximation has nothing to do with the value of the physical 
constants which are hxed by nature. It is solely decided on the basis of the particular form of the 
wave function. This in turn depends on the interaction potential and the boundary conditions of 
the problem. 

The purpose of this article is to derive the necessary and sufficient conditions on the inter¬ 
action potential and the boundary conditions under which the dynamical equations, namely the 
Schrodinger equation in the nonrelativistic case and Klein-Gordon equation in the relativistic case, 
are exactly solved by a semiclassical wave function. This is done in sections 2 and 3. Here the 
problem of the classihcation of all potentials which allow for exact semiclassical wave functions is 
solved. Section 4 includes a detailed analysis of the (1 -|- l)-dimensional Klein-Gordon equation. 
The results are then applied in section 5 for the study of solutions of the Wheeler-DeWitt equa¬ 
tion for FRW scalar field minisuperspace models. Here several exact semiclassical solutions are 
constructed. In section 6, the ideas and the results of the preceding sections are used to develop a 
novel semiclassical perturbation theory. The latter yields the semiclassical approximation in the 
zero-th order of the perturbation theory. The higher order corrections are shown to satisfy linear 
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differential equations with vanishing boundary conditions. In this way the information about the 
boundary conditions of the original problem is included in the zero-th (semiclassical) terms and 
the dehnition of the perturbation potential. The proposed semiclassical perturbation expansion 
is quite different from the traditional h and M~^ expansions used in quantum mechanics and 
quantum cosmology. 


2 Nonrelativistic QM: Schrodinger Equation 


Consider the Schrodinger equation 

, V'(O) = -00 (2) 

^ ^ [P - A(x; t)f + V(x;t) , 

where ^ is a state vector represented in the position representation by the complex scalar wave 
function (x|'^(t)) = t), A is an electromagnetic vector potential, and Id is a scalar interaction 

potential. 

Inserting Eq. ([^) in the Schrodinger equation @ and making use of (x|p = —zhV(x|, one 
obtains 


dtS{x;t) + H{x,p^;t)+ Q{x-,t) =0 , (3) 

9tp(x;t) + V • J(x;t) = 0 , (4) 

where H = if(x, p;t) is the classical Hamiltonian, p := B?, p* := VS", Q := —is 
the quantum potential, and J := pv* with v* ;= (p* — A) /m, is the probability current. 

Eq. (Q) is the quantum analog of the Hamilton-Jacobi equation 

dtS{x;t) + H{x,pp,t) = 0 , (5) 


of the classical mechanics, [|^. It differs from the latter because of the presence of the quantum 
potential Q. Eq. (H) is the continuity equation signifying the conservation of the probabilities. 

According to the above dehnition, the semiclassical or WKB approximation provides the exact 
solution of the Schrodinger equation, if and only if in addition to Eqs. (H) and (||), one has 




= 0 , 


( 6 ) 
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i.e., i? is a solution of the Laplace equation. In this case, Eq. reduces to the Hamilton- 
Jacobi equation (|^). Therefore, the necessary and sufficient conditions for the exactness of the 
semiclassical approximation are (^, (|), and (H). These equations may be viewed as three partial 
differential equations for the three unknown functions i?, S and V. 

Eq. does not involve time-derivatives. It is really a constraint equation which can be 
independently solved. Solving Eq. (§ and substituting the result in (^), one hnds a hrst order 
equation for v* which in turn yields p*. This leads to another hrst order differential equation for 
S. The potential V is then obtained by solving the latter equation and substituting the result in 
Eq. i). 

There is an alternative way of solving the continuity equation which involves writing it 
explicitly in terms of S, namely, considering the solution of 

V ■ (R^VS) = Rf , (7) 

where / := —2mdtR + 2VR- A-|-i?V ■ A. Now, let us dehne S := RS. Then, it is an easy exercise 
to show that S is the solution of the following Poisson equation 

\/^S = f. (8) 

Here I have used in addition to Eq. (0) the constraint equation (^. Hence, S is given by a 
solution of the Poisson equation (||) divided by a (non-zero) solution of the Laplace equation 
(1^). Note that both of these equations are second order elliptic differential equations with well- 
posed boundary value problems. Thus, R, S and consequently V are uniquely determined by the 
boundary conditions. This solves the problem of the classihcation of all nonrelativistic (scalar) 
quantum systems with an exact semiclassical solution for the Schrodinger equation by relating 
the latter to the boundary conditions of the Laplace and Poisson equations. It is also important 
to note that these boundary conditions may in general depend on time which appears in the 
corresponding equations as a parameter. 

In order to demonstrate the utility of these hndings in concrete terms, I shall next consider the 
case where the classical conhguration space is one-dimensional. Here one can pursue according to 
the former approach of integrating the continuity equation (^) by hrst solving for v*. 
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2.1 One-Dimensional Configuration Spaces 


Consider a quantum system whose configuration space is the interval [xi,X 2 ] C IR, and let the 
boundary conditions on the solution of the Schrodinger equation (|^) be given by t) = =: 

Ri{t) exp[iSi{t)/h], ip{x 2 ,t) = '02(^) =: R 2 it) exp[iS 2 {t)/h]. In this case, Eq. becomes d^R = 0 
which yields 


R = a(t)x + b(t) , 


(9) 


where 


_ Ri{t) - R 2 {t) _ xiRi{t) - X 2 R 2 {t) 


Xi -X2 Xi- X2 

Substituting (j^) in Eq. (^) and integrating the resulting differential equation, one obtains 


( 10 ) 


1 


n* = 


{ax + by 
S = d + m 


--cS,ax^-(cS,b+W,ay-2W,bx + c 


- {adtb + bdta)l 2 - 2bdtb h + c Jo 
o 


( 11 ) 

+ J Adx , (12) 


where c = c{t) and d = d{t) are functions of time determined by equating the right hand side of 


([T2|) with Si and S '2 at a; = xi and x = X 2 , respectively, and 

x^dx 


h : = 


{ax + by 


, fc = 0,l,2,3. 


More explicitly, one has 
In = 


a{ax + b) 


1 

b 

1 

b^ 

= - 

a^ 

In \ax + 0 H-- 

ax + b 

: -^2 Q 

ax + b — 2bm\ax + b\ H-- 

ax + b 


h = ^ 


1 b^ 

-{ax + by — 3b{ax + 6) + 3b^ In |aa: + 6| H- - 

2 ax + b 


The potential is then obtained using Eq. (H), namely 


m 


V = -dtS - -rvt . 


(13) 


It has the following general form: 

= C„(i) l„|a(«)x + K()] + - /a-4(^; 1)*. 

where Ci, with £ = 0,1, ■ ■ ■ 6, depend on a, b, c, and d. 

In view of the above analysis, one can reach the following conclusions: 
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— The condition of the exactness of the seniiclassical approximation together with the boundary 
conditions determine both the semiclassical wave function and the potential uniquely; 

— For xi —oo and X 2 —^ cxo, i.e., for a particle in IR, a smooth semiclassical wave function 
is not normalizable. It corresponds to a scattering state; 

— More general exact semiclassical wave functions may be obtained by allowing a countable 
number of discontinuities. The effect of these discontinuities is to divide the interval [xi,X 2 \ 
into a collection of subintervals in the interior of which the wave function and potential are 
given by the above expressions. The boundary conditions corresponding to each subinterval 
can be chosen freely. They determine the global structure of the wave function and the 
potential which can now be more complicated. This observation can also be used to devise 
an approximation scheme for the solution of the Schrodinger equation, by approximating 
the solution by a locally semiclassical one. 


Next, let us consider the following special cases: 

1) Constant Boundary Conditions: dt^Ji = dt'ilJ 2 = 0 with Ri ^ R 2 

In this case, a and b do not depend on time. This simplihes the above formulae considerably. 
One has: 


V 


c(i) 

(ax + by 


d(t) 


—dtd(t) + 


mdtc(t) 
a(ax + b) 


mc(t) 
a(ax + b) 
mc(ty 


2(ax + by 


+ J A(x)dx , 
J dtA(x)dx , 


where 


(14) 

(15) 


c 


li 


(S 2 - 


l(xi) , 


> 5 'i) + (71 - 72) 
Cl “ C2 


d = 


'y(x) := J A(x)dx , 


((1^2 — C2S1) + (C271 ~ C172) 

Cl - C2 

TTL 

Ci ■= ^ , with i = 1 , 2 . 

a(axi + b) 


In particular, for A = 0, 7 * vanish and c and d are constant. This leads to a further 
simplihcation of Eqs. O and ([^) and yields 


S 


d- 


mc 

a(ax + b) 


V = 


mvl 

~Y~ 


mC 

2(ax + by 
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In this case both the potential and the action function turn out to be time-independent. 
This corresponds to a semiclassical zero energy eigenfunction. 


2 ) Amplitude-Periodic Boundary Conditions, i.e., Ri(t) = R 2 (t) 

In this case, a = 0 and R = b{t). Then, a similar analysis leads to: 

-2bdtb X c 


n* = 


62 


S = d -\- m 


-{dtlnb) X + --r X 
0 


-\- Adx , 


V = —dfd — m 


— {dlhib)x‘^ + dt{c/b‘^)x 


m 


-26(9*6 X + c 
62 


— / dtA dx 


Here, c and d depend on the phases Si of -ipi according to: 


^ _ S 2 (t) - Si(t) ^ _ Si(t)(a 2 (t) - T.2{t)ai(t) 

a2{t) - ai{t) ’ Q; 2 (t) - Q;i(t) 

where 

■= Si -'fi + m{dt\nb) x^ . 

Note also that in this case, for A = constant, the potential V is at most a quadratic poly¬ 
nomial in X. For, b = Ri = R2 = constant, V is either a first order or a zero-th order 
polynomial in x. For example, for 6 = 1,5'* = Uit, with cn* being real constants, and A = 0, 
one has: 


c 

V 


(uJ2-C^l)t (UJ 1 X 2 - C02Xi)t 

m(x 2 — Xi) ’ X 2 — Xi ^ 

(- - - )[{^lX2 - ^2X1) + {UJ2 - UJi)x\ 

X2 - Xi 


(- )[{xiiX2 - UJ2X1) + {U2 - UJi)x\ 

X 2 - Xi 


1 

2 m 


(a ;2 - a;i)t 

{X2 - Xi) 


Another interesting case is when 6 = for some constant 00 . Then, the potential is 

a quadratic polynomial in x with the coefficient of the quadratic term being a constant, 
namely, —muj‘^/2. In the latter case if the phases Si and S 2 and the vector potential A are 
also time-independent, then so are c and d. Therefore, one has a time-independent quadratic 
potential, namely 

V = — — {ux + cY . 
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2.2 Multi-Dimensional Configuration Spaces 


For an n-dimensional configuration space, with n > 1, the classihcation of the exact semiclassical 
wave functions and the corresponding potentials is more involved. This is mainly because in this 
case the constraint equation is the n-dimensional Laplace equation = 0. 

In Cartesian coordinates, one can use the method of separation of variables to express i? as a 
sum of the basic solutions 


. (16) 

i=l 

where • • •, := x, Kj, a*, and are real functions of time which are determined by the 

boundary conditions on the wave function, Ki satisfy 


Y.’llfHit) = 0 , 


(17) 


Z=1 


with rji = ±1, and S (resp. C) stands for either of sin or sinh (resp. cos or cosh) depending on 
whether T]i = —1 or +1, respectively. Clearly, one can choose one of the r^j’s positive and others 
negative. 

Having found the expression for R, one then proceeds by integrating the continuity equation 
(^) which yields S. In view of the above analysis, S = S/R, where S' is a solution of the Poisson 
equation (|^). Using the well-known Green’s function methods of solving the Poisson equation, 
one may express S in an integral form. For example, if the conhguration space is IR^, then 




'3 


mdtR{x.'-, t) — Vi?(x'; t) ■ A(x'; t) — V ■ A(x'; t) 


X — X' 


8 ) 


where 5*0 is a solution of the Laplace equation determined by the boundary conditions. 

Similarly to the one-dimensional case, in the case that the configuration space is IR"^, a smooth 
semiclassical wave function cannot be normalized. More generally, it cannot vanish at inhnity, nor 
can it be localized. This is a direct consequence of Eq. (O- 

The main difference with the one-dimensional case is that here one has a much richer structure 
as far as the general form of the wave function and the potential is concerned. Unfortunately, 
since without knowing the specific form of the boundary conditions one cannot express R and S in 
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a closed form, an explicit classification of the semiclassical wave fnnctions and the corresponding 
potentials for n > 1 is not available. Nevertheless, it is evident that by choosing the bonndary 
conditions appropriately one can obtain a large variety of potentials. 

In order to demonstrate the validity of this claim, I shall next concentrate on the special cases 
where the amplitnde of the semiclassical wave fnnction is independent of x. A simple example of 
this is a particle in a cnbical cavity of side length L with the bonndary conditions: 

= N(t) , 

where the symbol d stands for the bonndary of the cavity. In this case, R = N{t) is the nniqne 
solution of the constraint eqnation (H), and Vi? = 0. This rednces Eq. (|^) to the simple Poisson 
eqnation 

V^S =-2mdt\nN , (19) 

where I have chosen the Conlomb gange so that V ■ A = 0. Note that the sonrce term on the right 
hand side of Eq. (pT]) only depends on time. Hence, one can dehne S := S + m[dtlnN] |xp and 
rednce this eqnation to the Laplace eqnation 

= 0 . ( 20 ) 

Since the set of solntions of the Laplace eqnation are in one to one correspondence with the set of 
bonndary conditions which is a very large fnnction space, one obtains a large class of potentials. 
Next, consider the following simple snbcases. 

1) Sq = Kit) ■ x|g, where /C is an x-independent vector-valned fnnction of time. Then, S = 
fC{t) ■ X clearly satishes the Laplace eqnation (p0|) and one has: 

S = /C(t) ■ X — In A(t)] |x|^ 

V = In A(t)]|xp — c}i/C(t) ■ X -|2m[c}tln A(t)]x + A(x;t) —/C(t)P . 

2m 

This case is the mnlti-dimensional analog of example 2 of section 2.1. 

2) Consider the case n = 2, i.e., a sqnare with bonndaries =: x = 0, L and x"^ =\ y = 0, L, 
and bonndary conditions: ^ = 0 for a; = 0, L, ?/ = 0, and S = a{t) sm.{nx/L) for y = L. 
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Then, one can easily show that the solution of the Laplace equation (^) is given by 


5 


a(t) sin(7rx/L) sinh^Tiy/L) 
sinh TT 


Hence, one has 


S 

V 


a(t) sin(7rx/L) sinh(7rj//L) 
sinh 71 

m[d‘f In A^(t)] |x 


— m[dt\n N{t)] 


2 9iQ;(t) sin(7rx/L) sinh(7rj//L) 


sinh TT 


2m sinh^ tt 


— cos(7ra:/L) sinh(7r|//L) — 2m[dt\n N {t)]x — 

L a{t) 


1 2 


+ 


— sin(7ra;/L) cosh(7r|//L) — 2m[dt\n N{t)]y — 

L a{t) 


n 2' 


These relations show that unlike the one-dimensional case, here the wave function and the 
potential can be quite complicated. 


So far, I have only considered cases for which Vi? = 0. This is precisely the condition de¬ 
manded in the traditional semiclassical approximation. It is usually argued that the semiclassical 
approximation is valid if the amplitude R of the wave function is a slowly varying function of 
X. As it is clear from the above analysis, this is only a sufficient condition, not a necessary one. 
There is an infinite number of examples where i? is a rapidly changing function of x but the 
semiclassical approximation is not only valid, but it yields the exact result. Specific examples can 
be constructed by simply choosing i? to be a rapidly changing solution of the Laplace equation. 
For instance consider a quantum system with the geometry of the preceding example, but the 
boundary conditions which lead to 


3) R{x,y;t) = R{x,y) = Nq sm{i7ix / L) smh.{i7ix / L) and S{x,y;t) = S{x,y) = 1C ■ x/i?, where 
Nq and £ are real and integer constants, respectively, and /C is a constant vector. One can 
easily check that R and S satisfy Eqs. (|) and (|^). However, 

Ittx, 




.Wtt 


COS 


,£'Ky 


l77X, 


,£'xy ^ 


(—)cosh(—)U/tan2 (—)+tanh ( —) 


can be made arbitrarily large by choosing large values for £, i.e., R is not a slowly varying 
function of x and y. Note also that in this case both R and S are time-independent. Thus, 
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the wave function describes a zero-energy eigenfunction of the time-independent potential 


^ +- IP - R -J ■ 

This example clearly shows how the present analysis generalizes the results of the traditional 
semiclassical approach to quantum mechanics. 


(yg)^ _ -1 

2m 2mi?2 


3 Relativistic QM: Klein-Gordon Equation 


Consider the Klein-Gordon equation 

10“ - A“)(a, - A,) - V(x)] i,(x) = 0 . (21) 

where are components of an electromagnetic gauge held, K is a scalar interaction potential 
(including the mass term in the massive case), and x stands for the four vector [x^). In the 
following, I shall follow the relativists’ convention for the Minkowski metric, namely, {rjnu) = 
diag(— 1 , 1 , ■ • •, 1 ), and set c = h = 1 . 

In the polar representation (|I|), the Klein-Gordon equation is written as 

- A>^){d^S - A,) + V + Q = 0 , (22) 

= 0 , (23) 

where Q := —d^dfj,R/R is the quantum potential and := p{d^S — A^), with p := is the 
conserved current. Again, Eq. ( [2^) is the quantum analog of the Hamilton-Jacobi equation for a 
classical relativistic particle 


{d^S-A^){d,S-A^)+V = Q, 


(24) 


and Eq. (|23|) is the continuity equation associated with the charge conservation. 
A semiclassical Klein-Gordon held is dehned by the condition 


Q: = 


d>^d.R 


df^d^R = 0 . 


(25) 


Therefore, the semiclassical or WKB approximation is exact, if and only if the relations (I 


(HI)) and (g^) are satished. As in the nonrelativistic case, these three equations may be used to 
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determine the three unknown real functions R, S, and V. This is done by hrst solving Eq. ([^) 
which is already decoupled from the other two. This is a wave equation for R. Its general solution 
is given by a linear combination of the functions of the form 

ly- = - X ■ k) (26) 


where x = {x^, x) belongs to the (n + l)-dimensional Minkowski space or a subset of 

and k is a unit u-vector dehning the null wave {n + l)-vector k = A;o(l;k)- A simple choice for 
which is essentially motivated by the Fourier analysis of the wave equation is the plane waves 
exp iik^x^ — X ■ k) . Once are chosen then the solution of Eq. (^5|) reduces to the determination 
of the coefficients of in the expansion of R. 


Next, one substitutes the expression for R in the continuity equation (^) and integrates the 
resulting equation. The basic strategy is similar to the nonrelativistic case. In terms of S the 


continuity equation (^) takes the form 


d^{R^d^,S) = RF , 


(27) 


where F : = 
dehning S 


2d^RA^ + Rd^A^. Eq. (^) is the analog of Eq. (|^. 


RS. This together with Eqs. (p5[) and (^) leads to 


It can further be simplihed by 


d^d^s = F , 


(28) 


i.e., S = S/R, where ^ is a solution of the inhomogeneous wave equation (^). Having found R 
and S', one can use the Hamilton-Jacobi equation (|^) to determine the form of the potential. 

There is a very important difference between the relativistic and nonrelativistic cases. Here the 
condition of the exactness of the semiclassical approximation leads to two second order hyperbolic 
equations, namely (|^) and (1^) , whereas in the nonrelativistic case one has two elliptic equations. 
One knows from the theory of hyperbolic differential equations that the boundary-value problem 
for such equations is not generally well-posed, i.e., for arbitrary boundary conditions, a solution 
may or may not exist and if it does, it may not be unique. The well-posed problem for a hyperbolic 
equation such as the wave equations (^) and ( p^) is the initial-value problem. In general for given 
initial data on a Cauchy hypersurface, one can solve these equations and determine i?. S', and V 
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uniquely. One of the consequences of the hyperbolicity of (|25|) is that unlike the nonrelativistic 


case, here R and therefore the semiclassical wave function can be localized. In particular, one can 
form a coherent wave packet which approximates the behavior of a classical particle. 

Restricting to the case where i? is a constant and adopting the Lorentz gauge = 0, one 

can reduce Eq. o to a (homogeneous) wave equation for S\ 

d^d^S = 0 , (29) 


The general solution of this equation is also given as a linear combination of functions of the form 


(1^) . This is sufficient to conclude that even for this special case S and consequently V can be quite 
complicated. This shows that there is a large class of potentials which allow exact semiclassical 
solutions of the Klein-Gordon equation. These potentials and the corresponding semiclassical 
Klein-Gordon helds depend in a crucial way on the boundary conditions^. This is especially 
important in quantum cosmology where there is an ongoing controversy regarding the choice of 
the boundary conditions for the wave function of the universe and also the form of the potential 
in the Wheeler-DeWitt equation. In particular, for the FRW scalar held minisuperspace models 
ii, the Wheeler-DeWitt equation is precisely a (1 -|- l)-dimensional Klein-Gordon equation of 
the form (pT]) . For these models the form of the potential is directly linked with the phenomenon of 
inhation BBi- On the other hand, most if not all the physical predictions which one hopes to 
derive from such models are relevant to the regions of the minisuperspace where the wave function 
is semiclassical. The results of this paper indicate that at least one can rule out the cases where 
the existence of a semiclassical solution (for some region of the minisuperspace) is inconsistent 
with the form of the potential (in that region). This together with the requirements imposed by 
inflation may be helpful in improving our understanding of quantum cosmology. This motivates 
a closer analysis of the Klein-Gordon equation in the (1 -|- l)-dimensional Minkowski space. 


^Here and in what follows, “boundary conditions” means “initial”, 


“boundary”, or “mixed initial-boundary 


conditions” for which there exists at least one solution. 
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4 Klein-Gordon Equation in (1 + l)-dimensions 


If the configuration space is the (l + l)-diniensional Minkowski space then the general solution 
of the wave equations ( ]25|) and (^) can be written in terms of four real-valued functions and 

III, 


R{x,t) = R^{x + t) + R_{x — t) , (30) 

S{x,t) = Sj^{x + t) + S_{x — t) + j dx'dt'G{xR]x'R')F{x\t') , 


where F is the same as the one appearing in Eq. (pTl) and G is the appropriate Green’s function 
for the one-dimensional wave equation, |^. The latter can be constructed out of the advanced 
and retarded Green’s functions given by G^{xR-^x'R') := [0{\x — x'\ ± (t — t')) — l]/2, where -|- 
and — label the advanced and retarded Green’s functions, respectively, and 6 is the step function, 
6{z) = 1, ii z > 0; 6{z) = 0, if z < 0. The usual choice in typical physical applications is 
the retarded Green’s function G~ which marks a particular direction of time. Note that if the 
electromagnetic potential is absent, F = 0, the last term in Eq. (|M|) drops, and there is no 
need for a Green’s function, in particular, a direction of time. Finally in view oi S = S/R one has 

1 


S{x, t) = 


S'+(a: -f t) -f- S'.(a: — t) -|- / dx'dt'G{x, t] x , t')F{x', t') 


. (31) 


F+ (x -1- f) -I- R- {x — t) [ 

Eqs. (|5I1|) and m together with Eq. (P^) show that in general the exact semiclassical wave 
functions and the corresponding potentials are classified by the set := {(F±, *S±)}, where C is 
the set of real-valued functions of a single real variable. 

Let us next concentrate on the case where there is no electromagnetic interaction. Then, in 
view of Eq. (^4[) the potential has the general form 

+ S^G^R' - S{R'^S' + R'G^)] 


V = 


(R, 


RA 


= -4 


f K 

K ] 

f S'_ 

R'_ \ 

[S+FS. 

R-f- -|- R— J 

U++^- 

F-|_ -|- R- J 


=; F[F±,S±] , (32) 


where R± = R±(x Ft), S± = S±(x F t), a prime means hrst derivative of the corresponding 
function. 




S+ + S_ 
R:+R_ ’ 


(33) 
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and the function(al) T is defined for futnre use. Eq. (^) is obtained by substituting (|30D and 


(|33|) in the Hamilton-Jacobi equation 

Next, consider the simple case R± = 1/2. Then, S = S^RS- and V = —4S^S'_. In particular, 
one has the following interesting examples: 


1) S is linear in time t: S± = uj±{x ± t)/2 for some constants uj± 

In this case, V = —a;+a;_ = constant. This includes the case of a free Klein-Gordon field of 
mass /i = a/— a;+a;_, since in this case V = jR. 

2) S is quadratic in t and x\ S± = h'±{x ± t)‘^/2 for some constants z/± 

In this case, one obtains a quadratic potential of the form V = — x‘^). This 

corresponds to a Klein-Gordon field with the time-dependent mass /i = 2^Vj^v_ t and 
quadratic interaction potential. The nonrelativistic limit of this case is a time-dependent 
harmonic oscillator with imaginary frequency. 

3) S' is a linear combination of exponential functions: S± = 

In this case, one has an exponential potential, V = where Vq := 

—Auj+uj-h'+h'-. Glearly, by choosing a;_ = ±a;+ =: ±a;, one obtains exponential potentials 
which depend only on t or x, namely V = and V = Voe^“*, respectively. 


These examples can also be described in the framework of the traditional semiclassical ap¬ 
proach, since R is chosen to be unity. Similarly to the nonrelativistic case, in order to demonstrate 
the generality of the present analysis, one must consider the cases where i? is a rapidly changing 
function of t and x, but the semiclassical approximation is nevertheless exact. Again typical ex¬ 
amples can be constructed starting from a rapidly changing solution of the wave equation (^5[) . 
For instance, consider the case 

4) R = and S = 

Then one can check that Eqs. ( pS]) and (pT]) are satisfied and the corresponding semiclas¬ 
sical wave function is an exact solution of the Klein-Gordon equation defined by the time- 
independent potential V = — 
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Eq. (|32D provides a classification of the potentials for which an exact semiclassical solution of 
the (1 + l)-diniensional Klein-Gordon equation, with = 0, exists. In practice, however, it is the 
potential which is given, not the wave function. Hence, a more interesting question is whether for 
a given potential V = V{x, t) there is a set of boundary conditions which yields an exact solution 
of the Klein-Gordon equation. One can alternatively ask whether a given potential belongs to the 
image of the function JF dehned in Eq. (|52|). 

In order to answer these questions, I shall hrst consider a special class of boundary conditions, 
namely = 0. For this class one can show that T is not onto and the set of potentials of the 


from Eq. (|3^) , with = 0, forms a small subset of all possible potentials. To see this, let us hrst 
substitute = 0 in (|3^) . The resulting equation may then be viewed as a differential equation 
for 5_, while and 5+ are treated as undetermined functions. This leads to 


p/ _ p2 




4 / 2'^4 


= 0 , 


(34) 


which is a consistent hrst order ordinary linear diherential equation for S- provided that the 
bracket on its left-hand side depends only on x — t. This puts an strong restriction on the form of 
the allowed potentials. Namely, the potential must be of the form: 


V (x, t) = X{x — t)S'j^{x -|- f) -|- y{x — t)Sj^{x -|- t)S'j^{x 1) , 


(36) 


where X is an arbitrary real-valued function and y := 4i?'_/i?T 


Proposition 1: Let u := x — t and v := x-\-t be null coordinates in and V : ^ JR be 

an analytic function at (0, 0) G Then, a necessary condition for V to satisfy Eq. 
for some functions T : IR —IR, 3^ : IR —IR, and S'+ : IR —IR, is that the coefficients 
Vjn in the power series expansion of V must satisfy one of the following two 

relations 




^jmi 


Vkim2^j'in2^^k2n i\^k2m2^jm2 ^^ 2^1 




^^im2 ^k2mi 


^jmi '^kin 


(36) 

(37) 


with (j, ^ 1 ,^ 2 ) and (n, 777 . 1 , 777 . 2 ) being triplets of mutually diherent arbitrary non-negative 
integers. 
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Proof: Substitute the power series expansions 


P = Viu,v) =: ^ P,- 


jn 




WV 


j,n=0 


X = X(u) =: 

j=i) 


y = J>(«) =: , S+ = S+(!.) =: Y.Sn 


(38) 

(39) 




n=0 


in Eq. 


This leads to 


n—1 


VjYi ^ ^ ^“1“ mi ("^“1“ 1) ^ ^ ('H/ T/l) *5*7^2_|_i5*72—m : 

m=0 m=0 

(40) 


where .Ej := Xj + S'o3^j. Next, solve for the sum on the right-hand side of (|40|). The result is 


n—1 


^ ^ Sm+lSn—m {^jn T ^') ■ 


m=0 


Clearly the left-hand side of this equation is independent of j. Hence, its right-hand side 
must also be independent of j. Writing the right-hand side for two different values of j and 
equating the results, one has 


c -1 

" n[y, yk. 


-1 


1 hfc,n—1 




(41) 


,3^. yk J 

Next let US express Sn using two different values of k, say ki and ^ 2 . Equating the two 
expressions and simplifying the result, one obtains 

yk2yjn yjyk2n 


y - 


+ (yki^k^ — yk^^ki] 


(42) 


Since n does not appear in this equation except in the square bracket, the content of the 
square bracket must be independent of n. This argument may be used to determine by 
equating the square bracket on the right hand side of with its value for n = mi. This 
leads to 

3^. = 3^A:2W,n , (43) 


where 


W,. : = 


(cVfc2mi ^^imi)^'n (cVfc2n 


jmi 




kimi ^k2n 


(44) 


and c := yki/yk 2 - Once again, Wjn = yj/yk 2 must be independent of n, i.e., for all m and 
'n, yVjn = y^jm- In particular Wjn = Wjm 2 , where m 2 is some arbitrarily chosen fixed non¬ 
negative integer. This equation may be used to express in terms of , Vjm 2 , I4in, Vk^n, 
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and c namely 




cV, 


k2n 


V, 


kin 


cVu 


k2m2 


K 


kim2 


^k2mi ^kin 


^kimi ^k2n 


cV^2mi ^kimi \CV^2mi ^kimi 
'^kimi^k2n 


^ ^^2^1 ^kim2 


^kimi ^^2^2 , 


'^jmi “t“ 




jm2 • 


Next, consider the following two possibilities: 


(45) 


I) C ^ l^imi/ldc2 mi • 

Then, the right-hand side of ( ^5D does not actually depend on c. In this case, one finds 

{y^kimi^jmi ^kim2^jm2)^k2n T (h^2m2^"m2 ^k2mi^jmi^^kin 


Vjn = 


h^imi h/i:2^T?,2 ITim,oId 


kim2 ^k2mi 


which is just Eq. 
ff) C ^di;imi/h^2 mi • 

Then, according to the definition of c, i.e., c := yki/yk 2 and Eq. (0) either all Sn 
vanish — this corresponds to the trivial case Vjn = 0 — or Xj/yj = X^jy^^ i.e., the 
ratio Xj/yj does not depend on j. The latter implies X = rjy for some constant r]. 
Moreover, in this case Vjn must satisfy 

T/ _ yjmiVkin 

^jn — .. , 

r kimi 


which is just Eq. (0)- ° 


Note that in the latter case, Eq. (^7|) together with X 


rjy and Eqs. ( 140|) and (^) lead to 


V,n 

V{x,t) 


n 

Xnyj y Xn ■ Tji^n T IjiSfi-i-i T ^ : (n TTI T X'jSniSn—m y 

m=0 


[t] + S+{x + t)]S'^{x + t)y{x - t) . 


(46) 

(47) 


Now, substituting Eq. (^) in Eq. (|3^ , one can easily show that indeed 5*+ drops out of this 
equation and one obtains S- = r]. Therefore the exact semiclassical wave function ^jJ = is 
given by i? = and S = {rj + S+)/R^. For example, consider choosing 

5) r = and (r/ + 3^)3'^ = for some constants fj.± and oj±. 

Then, one has V = In particular for a;+ = —a;_ =: a;, the potential 
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depends only on t. In this case, one has 


1 / = , 


R = 




I V. 


2uj 


S = ±2^ 


+ v+e 


—LOt 


11-6 


+ i/_e 


— Ljt ’ 


(48) 


where z/± are also constants. The appearance of the square roots in these equations is an 
indication that for certain choices of /r± and either Rot S can become imaginary in some 
regions of the Minkowski space. Since R and S are assumed to be real, such a semiclassical 
solution does not exist in these regions. 


This is another example of a case where the amplitude of an exact semiclassical solution is not a 
slowly varying function of its arguments. 

A simple consequence of Proposition 1 is: 


Corollary: The potentials of the form which are analytic at (x = 0,t = 0) form a 
proper subset of the set of all potentials, i.e., for an arbitrary potential which is analytic 
at (x = 0,t = 0), an exact semiclassical solution of the (1 + 1)-dimensional Klein-Gordon 
equation with R+ = 0 and = 0, may not exist. 


Eqs. (1^ and (|^) provide a useful criterion for Ending out whether a given potential allows 
for an exact semiclassical solution with = 0 or not. If the result is positive, then Eqs. O. 
(|4^), and (^31), may be used to determine and S+. These are then used to integrate Eq. 


which yields S- and consequently the wave function -0 = exp[z(S'_ + S'+)/i?_]. 

Proposition 1 only applies to the cases where R+ = 0. One might try to employ a similar 
method to treat the more general case, where is also an undetermined function. For this 
purpose, one must first write Eq. ([3^ ) as a polynomial equation in V, R±, S± and their derivatives 
and substitute the power series expansions of these functions in the resulting expression. This leads 
to an infinite system of very complicated coupled nonlinear algebraic equations for the coefficients 
of R±, S± whose analytic solution has not been possible. Although a similar proof is lacking for 
the most general case, further inspection of Eq. (|3^ suggests that this equation also restricts the 
form of the potential. Hence, in general for an arbitrary potential an exact semiclassical solution 
of the Klein-Gordon equation does not exist. 
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5 Exact Semiclassical Wave Functions of the Universe 


Consider the Wheeler-DeWitt equation for the closed FRW cosmological model coupled to a scalar 
field 0, 0, |], 


-dl + al + 


e6“V(0) ^/' = 0, 


(49) 


where a := In a, a is the scale factor of the FRW model, V is the matter held potential. Here, 
the cosmological constant is assumed to vanish, the usual factor ordering prescription [0, Q is 
adopted, and the natural units in which the Planck mass is set to unity is used. Clearly, this is a 
Klein-Gordon equation in (1 + l)-dimensional Minkowski (minisuper)space with potential 


V = -e^“ + e®“V(0) , 


(50) 


where (a,0) play the role of {t,x). 

The solutions of Eq. (^) have been studied mostly for the massless (V = 0) and massive 

0, |T^ H, This is done by making use of 


(V = nrcjr) scalar helds in the literature |]^ 


different approximation schemes except for the rather trivial and much less interesting massless 

are 


case for which the exact solution is known [0^, The approximate solutions of Eq. 
usually developed by making particular assumptions for the boundary conditions, semiclassicality 
of the solution, or restricting to particular regions of the minisuperspace in which the Wheeler- 
DeWitt equation (|49|) simplihes. 

In view of the developments reported in the preceding sections, the assumption of the exactness 
of the semiclassical approximation provides a direct link between the choice of the boundary 
conditions and the form of the potential. This is done through the function JF dehned by ([3^ ) 
which can be viewed as a function from the set of boundary conditions to the set of potentials 
V. In section 4, I have shown for the case R+ = 0 that this function is not onto, i.e., there are 
potentials which do not admit exact semiclassical solutions. Here, however, one is interested in a 
special class of potentials, namely those of the form (|50D . Hence, the relevant problem is to hnd 
the intersection U of the image of JF and the set of potentials of the form (^). One can easily 
show that U is not empty. For example the massless case, where V is an exponential function of 
a, can be easily put in the form (^), i.e., it belongs to U. In fact, two possible choices for R± and 
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S± which lead to this potential were already given in examples 3 and 5 of section 4 (with a; = 2 ). 


The following are nontrivial examples of the potentials of the form (|50|) which also belong to U. 
They are obtained by setting = 0 and making simple choices for the functions X, y, and 5*+ 
of Eq. so that Eq. is also satished. This together with Eq. yield S-. 


1) V = with A > 0 

In this case, the choices X = — R_ = 2e'^““/\/A, and satisfy 

Eq. (pSl) . The solution of (|3^) then leads to + 1/A, where c is a constant. The 

exact semiclassical solution is given by 


R = 


2e^- 


Vx 


, ^ = yA(e'^+^“ + 


2 e° 


A 


+ 2 c) ; 


2 ) V = Ac with arbitrary A 

In this case, one has X = R- = 2 ci[c 2 — 5_|_ = Cie‘^+“, and 

Aci 


= 


2g-3(<A-a) ^ 3g2e-('^-“) + 3c^e‘^- 


tan ^ J — 1 + cs 


/c^^2{4>-a) _ X 

where Ci, C 2 and C 3 are constants. The exact semiclassical solution is therefore given by 

2ci 


R = 


S = 


C2 - 


3</’ + “ A 


2e-3(<>-“) + 3c2e-('^-“) 


\lC2 — ^ ^C2e'^^X>-o‘) — 1 + 


16 


C 3 


Note that in this case, there is always a region of the minisuperspace dehned by “ < C 2 
where a semiclassical solution does not exist. 

The next logical step is to explore the existence of exact semiclassical solutions of the Wheeler- 
DeWitt equation with matter potentials of the form V = A(/)^^. These are among the potentials 


which lead to inflationary classical solutions, |]^. The simplest case is that of a massive scalar 
field, i.e., A = nR, p = 1. In the remainder of this section, I shall restrict to the semiclassical 
solutions with = 0. The existence of this type of solutions can be easily decided using Eqs. (|36|) 
and (|37|). One simply needs to compute the coefficients Vjn of ( PB[ ) and check whether they satisfy 
one of these equations. A simple calculation shows that for p = 1, 2 , 3, none of these equations 
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are satisfied. Hence, the matter potentials A0^, A0^, and A0® do not admit exact ‘right-going’ 
(/?+ = 0) semiclassical solutions. This is done by choosing the integers j, n, ki , k 2 ,mi, and m 2 
in such a way that both Eqs. (^) and (^Tl) fail. In order to demonstrate this, let us denote the 
right-hand sides of Eqs. (|36[) and (^) by and , respectively. Then, 


for p = 1 and ki = mi =0, ^2 = m 2 = 1, one hnds 


1/ - 4 + - 4 + - 

r 22 — — 4 + — , \/22 — — 4 H -^-— , ^22 — 




32 3 + 4 , 32 3 ^2 16 ’ 3 


llA 3A2 


for p = 2 and ki = mi =0, ^2 = m 2 = 2, one hnds 


f"31 = 


8 A 


y(l) _ T/(2) _ _ 

^ ) *"31 — *^31 — 3 ) 


33 g g , 33 9 + 8 ’ 9 


for p = 3 and ki = mi = 0, ^2 = m 2 = 3, one hnds 


f"24 = -V 


15A 


y(l) _ T/(2) _ 
>^24 — ^24 — 


4 

’3 


Note that for A = 0, i.e., the massless case, the values of Vjn in the above list match the values of 
and Vj^\ This is in agreement with the fact that for the massless case one does in fact have 
right-going exact semiclassical solutions. Eqs. (H) with (t,x) —>• (a,0), a; = 2, and p_ = —l/p+, 
provide a concrete example. 

Furthermore, one knows from the studies of the inhationary cosmological models that for the 
polynomial matter potentials A0^^ the coupling constant A must be a very small number. For 


example for the massive case, where p = 1 and A = m^, these theories predict m ~ 10 


-6 


i.e., 


A ~ 10 [13]. Thus, although there is no right-going semiclassical solutions for p = 1,2,3, at 


least for small values oi u = (p — a and v = (p + a, where one can neglect forth and higher order 
terms in the power series expansion of V, the semiclassical approximation seems to be reliable. 
The phrase semiclassical approximation is used in quantum cosmology in a very crude way. One 

in to reduce 


usually makes additional assumptions such as the adiabaticity of the evolution 
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the situation to the one-dimensional quantum mechanical case. In this way, one is able to express 
the condition of the validity of the semiclassical approximation as a simple limitation on the range 
of values of the matter potential, 0. The consistency of these assumptions with the validity of 
the semiclassical approximation is either left unchecked or a set of conditions is imposed which 
render the scheme consistent. These conditions are usually sufficient conditions not necessary. 
Hence, in general they may be too restrictive. The situation is very similar to restricting the exact 
semiclassical wave functions to those with slowly varying amplitudes. As shown in the preceding 
sections, this is an absolutely unnecessary restriction. The approach pursued in the article also 
allows for a precise definition of a more general semiclassical approximation where the solution to 
the dynamical equations is approximated with the general semiclassical wave functions introduced 
in section 1. This will be discussed next. 

6 Semiclassical Perturbation Theory 

Let us hrst define a semiclassical potential Vq to be a potential which corresponds to an exact 
semiclassical solution of the dynamical equation. In view of the results of sections 2 and 3, the 
set of semiclassical potentials is in fact much larger than one usually expects. This suggests a 
generalized notion of semiclassical expansion and in particular semiclassical approximation which 
corresponds to a perturbation theory around the semiclassical potentials. 

Let V be an arbitrary potential, ip = be the solution of the dynamical equation, and 

e G IR be a perturbation parameter. Then, the semiclassical perturbation theory corresponds to 


V = 

Vo + 6V , 

6V = eVp 

(51) 

R = 

Rq 6R , 

OO 

5R = Y.e^R, , 

(52) 



£=i 


S = 

So + 5S , 

OQ 

5S = Y.e^St, 

(53) 


£=1 


where Vq and ipo = are the semiclassical potential and wave function defined by the 

boundary conditions of the problem. Substituting Eqs. (|^) - in the dynamical equations 
and treating e as an independent parameter, one obtains an infinite family of equations which can 
be iteratively solved to yield Ri and S^. In particular, the equations obtained in the £-th order are 
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two linear coupled differential equations in Ri and with vanishing boundary conditions. This 
is because the original boundary conditions are already imposed in the determination of Rq and 
5o. 

Eqs. and (|5^ yield what one might call a semiclassical perturbation expansion. Note 


that they are not power series in h but in the perturbation parameter e. The zero-th order terms 
correspond to the semi-classical wave function. Thus, the semiclassical approximation is dehned 
by i? ~ i?o and S ^ Sq. Note also that the semiclassical wave function 'i/jq and potential Vq are 
uniquely determined by the boundary conditions. The choice of the perturbation parameter is, 
however, made by the physics of the problem. A typical example of a perturbation parameter is 
the coupling constant A of the preceding section. 

Let us next list the equations governing the hrst and second order terms in the semiclassical 
perturbation expansion. 

— Nonrelativistic QM: Schrodinger Equation 

First order (post-semiclassical) corrections, i.e., equations determining Ri and Sp. 

dtSi + -{VSo-A)-VSr + V^-^^^^ = 0 , 

m zmRo 

1 


2dt{RoRi) +-V ■ Rt,VSi + 2RoRi{VSo- A) 
m L 


= 0 . 


Second order corrections, i.e., equations determining i ?2 and S' 2 : 

1 r__ _ .. _„i 1 fV^R 2 RiV^Ri 


dtS 2 + — {VS,y + 2(VS'o - A) ■ V^2 

2m L 


2m \ Rn 


Rl 


2dt{RoR2) + dtRj + -V ■ (Rl + 2 RoR2){VSo - A) + 2RoRiVSi + RIVS 2 

m L 


= 0 , 

= 0 . 


These equations are obtained by substituting Eqs. (0 - (|53D in Eqs. (|) and 


Relativistic QM: Klein-Gordon Equation {c = fi = 1) 


First order (post-semiclassical) corrections, i.e., equations determining Ri and Sp. 


2(a^S'o - A^d^Si + Ep + 


Rn 


d„ 


2RoRyd>^So - Ay + Rld>^Si 


= 0 , 


= 0 . 


(54) 

(55) 
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Second order corrections, i.e., equations determining i ?2 and S' 2 : 


2(a^S'o - A^^)d^S2 + + 




Rn 


Rl 


d„ 


Rld^S 2 + 2 R^Rid^^Si + + 2R0R2WS0 - 


0, 

0 . 


These equations are obtained by substituting Eqs. (|^ ) - (|53D in Eqs. (P2|) and (P^l). They 
can be further simplified. For example, consider the first order equations ( |5^ and (^3]), and 
define pg •= — A^ and T±i := Ri± RqSi. Then using the fact that Rq and Sq define a 

semiclassical wave function, i.e., d^d^Ro = 0 and 9^(i?oPo) = 0) can show that Eqs. 


and (|55[) are equivalent to 


± 2p^{d^ - In i?o)] T±i = -RoV, 


p • 


(56) 


These are two separate equations for T±i whose solution yields Ri and Si. 

These equations appear to be even more difficult to solve than the original dynamical equations. 
Note, however, that they are to be solved with vanishing boundary conditions. The advantage of 
this scheme is that the information on the boundary conditions of the original dynamical equation 
is restored in the definition of Rq and Sq. The higher order corrections are affected by these 
boundary conditions only through Ro, So and the perturbation potential Vp which appears in the 
first order of perturbation. In this sense, the semiclassical perturbation theory has a universal 
character. 

As a concrete example consider the Wheeler-DeWitt equation of section 5 with a matter po¬ 
tential of the form V = A/(0), where / is some real-valued function. Let the boundary conditions 
be such that one recovers example 5 of section 4 with uj = 2, ii± = =pl, z/_ = 0, and = v. Then 
replacing (t, x) with (a, 0), (V, i?, S) with (Vg, i?g, S'g), and adopting the positive sign in Eqs. (^8[), 
one has 

Vo = -e^“ , Ro = , ^g = 2^. (57) 

These correspond to an approximate semiclassical solution which is an exact solution of the mass¬ 
less case. Clearly, the perturbation potential Vp is given by e®“/(0) and the perturbation param¬ 
eter is A. The first order (post-semiclassical) correction (i?i. Si) to (i?g, S'g) is obtained by solving 
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Eqs. (|56|) which take the form 

{-So + ± [«So - 1) + m + 1)1} r±i = . (58) 

with 

- 2e-*) 

" ^ ^^g2(^-Q) _ g4^. 

Although solving Eqs. (|58|) seems much more difficult than solving the original Wheeler-DeWitt 
equation (^), one must recall that these equations are to be solved with vanishing boundary 
conditions. In this way one can at least devise an efficient numerical scheme which can treat the 
solution of these and the equations for the higher order corrections for arbitrary matter potentials. 

In general, such a scheme should hrst compute the semiclassical potential Vq and wave function 
Roe‘^^° using the boundary conditions. This would yield the perturbation potential Vp (after the 
perturbation parameter is identified according to the physical characteristics of the problem). 
Then, it should numerically integrate the equations satisfied by Ri and Si with vanishing boundary 
conditions. 

Finally, let me emphasize that for the Wheeler-DeWitt equation with the massive scalar held. 
The perturbation parameter is already extremely small A = = 0(10“^^), therefore the hrst 

order corrections provide solutions which are valid up to the order = 0(10“^^). This suggests 
that the domain of the validity of the hrst order perturbation theory is indeed quite large. 

7 Conclusion 

In this article I have tried to demonstrate how the simple observation that the traditional condition 
of validity of the semiclassical approximation is only a sufficient condition, can be used to introduce 
the notions of exact semiclassical wave function and potential. 

I have shown that the semiclassical wave function and potential are both determined by the 
boundary conditions of the problem uniquely. I have also given the full classihcation of exact 
semiclassical wave functions and potentials for the Schrodinger and Klein-Gordon equations in 
arbitrary dimensions. The analysis of the one-dimensional Schrodinger equation is much simpler 
than the multi-dimensional case. For the Klein-Gordon equation in (l-l-l)-dimensions which is 
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directly relevant to the solution of the Wheeler-DeWitt equation for FRW scalar field minisuper- 
space models, I have explicitly constructed semiclassical wave functions and potentials. For this 
case, I have also developed a practical criterion for checking whether a given potential allows for a 
right-going exact semiclassical solution. I have then used this criterion to study the semiclassical 
solutions of the minisuperspace Wheeler-DeWitt equation. For the polynomial matter potentials 
of the form with p = 1,2,3, I have shown that a semiclassical solution of the Wheeler- 
DeWitt equation does not exist. However, the non-existence proof relies on the non-zero value of 
the coupling constant A which is expected to be an extremely small number. This motivated the 
development of a semiclassical perturbation theory which yields the semiclassical approximation 
as the zero-th order term in the semiclassical perturbation expansion. The higher order terms 
satisfy coupled linear differential equations with vanishing boundary conditions. 

The resulting semiclassical approximation and the domain of its reliability are different from 
the traditional semiclassical approximation. They coincide for the cases where the semiclassical 
approximation is exact, i.e., the potential and wave function are semiclassical. The semiclassical 
expansion developed in this paper is a perturbation expansion about the exact semiclassical po¬ 
tential dehned by the boundary conditions. This is in contrast with the traditional semiclassical 
expansion which is an h or ^-expansion of the solution of the dynamical equation. 
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